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ABSTRACT 

We calculate the radiation spectra from advection-dominated accretion flows (ADAFs), taking into 
account the effects of a global magnetic field. Calculation is based on the analytic model for magnetized 
ADAFs proposed by Kaburaki, where a large-scale magnetic field controls the accretion process. Adjust- 
ing a few parameters, we find that our model can well reproduce the observed spectrum of Sagittarius 
A*. The result is discussed in comparison with those of well-known ADAF models, where the turbulent 
viscosity controls the accretion process. 

Subject headings: accretion, accretion disks — radiation mechanisms: non-thermal — Galaxy: 
center — magnetic fields 



1. INTRODUCTION 

The optically thin, advection-dominated accretion flows 
(ADAFs) have been studied by a number of authors dur- 
ing past several years (e.g. Narayan & Yi 1994, 1995a, 
b; Abramowicz et al. 1995; Nakamura et al. 1997, Man- 
moto, Mineshige & Kusunose 1997; Narayan et al. 1998). 
These models are very successful in describing both spec- 
tra and dynamics of accreting black hole systems such as 
those in binaries and in low-luminosity active galactic nu- 
clei (AGNs). The observed spectra can be explained as 
follows. The radio emission is due to the synchrotron 
emission in turbulent magnetic fields in the accretion flow. 
These synchrotron photons serve as seed photons for the 
inverse Compton process by hot electrons. Once-scattered 
Compton photons are mainly distributed in the optical 
band and twice-scattered Compton photons, in soft X- 
ray band. Bremsstrahlung due to electron-electron and 
electron-proton collisions gives rise to the observed hard 
X-ray spectra. 

Thus, these ADAF models provide a good framework 
for understanding the observed spectra. In these models, 
both angular momentum transfer and energy dissipation 
in the accretion flow is assumed to be undertaken by the 
turbulent viscosity whose size is specified by so-called a 
parameter. For this reason, hereafter we call this type of 
models the "viscous" ADAF model in this paper. The 
magnetic fields are regarded as of turbulence origin and 
are described by another parameter [3 which specifies the 
ratio of the magnetic pressure to the gas pressure. 

However there is no reason to believe that the turbu- 
lent viscosity is the only candidate that controls the ac- 
cretion processes. Rather, it is quite natural to think that 
some types of global magnetic fields may play an essential 
role. Indeed, there are some evidences for the presence 
of such an ordered magnetic field in the central region of 
our Galaxy (e.g., Yusef-Zadeh, Morris & Chance 1984). 



As Kato, Fukue & Mineshige (1998) has pointed out the 
hydromagnetic turbulence in accretion disks may also gen- 
erate global magnetic fields by dynamo processes due to 
the presence of helical motions. In view of such circum- 
stances, another type of ADAF model has been proposed 
by one of the present authors (Kaburaki 1999, 2000; here- 
after referred to as K99 and KOO). In order to distinguish 
it from the above viscous ADAF models, hereafter we call 
it the "resistive" ADAF models since energy dissipation 
in the accretion flow is due to the electric resistivity and 
angular momentum transfer is supported not by the vis- 
cosity but by the magnetic stress of a large scale magnetic 
field. 

The purpose of the present study is to calculate the 
expected radiation spectra from ADAFs in a global mag- 
netic field based on the resistive ADAF model, in order to 
compare its predictions with those of the viscous ADAF 
models. As a most suitable candidate for such a compari- 
son, Sgr A* is taken up here because it has been observed 
in many wave lengths as the nearest galactic nucleus and 
its spectrum has been reproduced many times by the suc- 
cessively advancing viscous ADAF models. 

In §|[ we introduce the set of analytic solutions for re- 
sistive ADAFs in a suitably scaled form and discuss their 
basic characteristics. The relevant radiation mechanisms 
and_the methods of calculation of the fluxes are jdescribed 
in §kI These schemes are applied to Sgr A* in §W and the 
resuTts are discussed in comparison with those of the vis- 
cous ADAF models. Finally in §pl we summerize the main 
results and discuss some related issues. 



2. RESISTIVE ADAF SOLUTION 

As a basis of our calculation of spectra, we introduce 
here the set of analytic solutions constructing the resistive 
ADAF model. This set may be considered as a counter- 
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part of that found by Narayan & Yi (1994, 1995a) in the 
viscous ADAF scheme, but it should be emphasized that 
the former is not a self-similar solution as the latter. In 
the resistive ADAF model, there arc three basic quantities 
and one parameter: mass of the central black hole M , mass 

accretion rate M, strength of the external magnetic field 
I Bo I and half-opening angle of the flow A, respectively. We 
introduce the following normalizations for these quantities 
and for the radial distance R: 

M . M , |Bo| . A R 

(1) 

where i?out denotes the radius of the disk's outer edge. The 
Eddington accretion rate is defined by Me = L-^/ (Q.lc^), 
which includes the efficiency factor of 0.1. Note that this 
definition of M-^ and the normalization factor for black 
hole mass are different from those in KOO. The latter is 
chosen here as 10^ Mq for the convenience of the discus- 
sion of Sgr A*. 

In spherical polar coordinates, the radius of outer edge 
and the radial-part functions of relevant physical quanti- 
ties (for their angular parts, see also K99, KOO) are written 
as 

i?out = 1.5 X 10^^ ho'^'^m^'^m^/^ cm, (2) 
\b^{r)\ = lQ5-hor-^ G, (3) 
WK(r) = 0.9 X 10*^ bl'^^fn-^'^m^'^r-^/^ cm s^\ (4) 
P{r) = A.Q 5-%lr-'^ dyne cm" ^ (5) 
T(r) = 1.8 X 10Uo^'m-2/W/-V-i K, (6) 

p{r) = 1.3 X 10-15 5-%y^m"''m-^"'r-^ g cm'^, (7) 

where is the toroidal magnetic field, Uk is the Kepler 
velocity, P is the gas pressure, T is the temperature (com- 
mon to electrons and ions) and p is the density. Owing 
to the non-negligible pressure term in the radial force bal- 
ance, the toroidal velocity in the disk is reduced by a factor 
of l/^/2> from the Kepler value. The surface density and 
the optical depth are given, respectively, by 

E(r) = E = 4.1 S-^'^tn^'^m^'^ g cm'^, (8) 

Tes(r) ~ ^/^esS = 8.2 X lO'^ , (9) 

where the opacity for electron scattering. Note that 
these are independent of r in the present model. 

Fig. 1 shows a schematic picture of an accretion disk 
in a global magnetic; field, whose precise structure is de- 
scribed by the solution given above. Otherwise uniform 
external magnetic field is twisted by the rotational motion 
of accreting plasma, and there develops a large toroidal 
magnetic field in the middle latitude region. The behav- 
ior of this component, especially within the geometrically 
thin accretion flow, is given as — 6,^tanh^ in the resistive 
ADAF solution, where ^ = (6' — 7r/2)/A is the normalized 
angular variable. Owing to the appearance of a global 6^, 
angular momentum of the accreting plasma becomes able 
to be carried away by the magnetic stress to distant re- 
gions along the poloidal magnetic field. The extraction 
of angular momentum guarantees the inward motion of 
the plasma, which gradually becomes large until it reaches 
near the rotational velocity at around the inner edge of the 
accretion disk. Although the magnetic lines of force are 



also bent inwardly toward the center of gravitational at- 
traction, it has been shown that the dominant component 
is the toroidal one. This component also plays an essen- 
tial role in the plasma confinement toward the equatorial 

plane and keeps it geometrically thin through its magnetic 
pressure. 

Before going into the detailed discussion of the radia- 
tive processes, we briefly mention some similarities and 
differences in the basic features of the viscous and resistive 
ADAF models. It is worth noting that, in spite of the es- 
sential difference in the mechanisms of angular momentum 
transport and energy dissipation, the predictions for quan- 
tities such as temperature, density and optical depth are 
quite similar in both models. The temperatures in both 
models are as high as a fraction of the virial temperature of 
ions. Indeed, this is the case in the viscous models, though 
the electron temperature may deviates from it in the in- 
ner regions (e.g., Narayan & Yi 1995b), and so also in the 
resistive model as can be confirmed from the analytic ex- 
pression of T (K99, KOO). Such a high temperature makes 
a sharp contrast with the case of standard a-disks (e.g., 
Frank, King & Raine 1992). Further, for a sub-Eddington 
mass accretion rate, the optical depth is dominated by the 
electron scattering and is smaller than unity. These are 
therefore common features of sub-Eddington ADAF mod- 
els as expected. 

As one of the main differences, it may be stressed that 
the magnetic field in the resistive ADAF model is an or- 
dered magnetic field and is determined self-consistently in 
the model from a boundary value. Therefore, the strength 
of the field is not a parameter as in the viscous ADAF 
models. The ordered magnetic filed extract angular mo- 
mentum from the accreting plasma and confines it in a disk 
structure against the gas pressure. Gravitational energy is 
released in the disk as the Joule heating and also as com- 
pressional heating of the fiow. In the above analytic model 
of resistive ADAFs, these energies are fully advected down 
the stream. We calculate the radiation from the disk as a 
small perturbation from this solution. 

Another distinction may be in the energy partition be- 
tween the electron and ion components. The viscous 
ADAF models assume that the viscous dissipation, which 
is large at large radii, heats mainly ions. Since the effi- 
ciency of radiation cooling is very small for ions compared 
with electrons and since energy transfer to the electron 
component is estimated to be negligible (Manmoto et al. 
1997) except in outer portions, the flow is fully advective in 
most portions. For the electron component, on the other 
hand, the radiative cooling is balanced by the advective 
heating, near the inner edge. The electron temperature, 
therefore, deviates downwards largely from the ion tem- 
perature thus realizing a two-temperature structure. 

In contrast to the viscoTis heating, the resistive dissi- 
pation becomes large at small radii and seems to pref- 
erentially heat the electron component as suggested by 
Bisnovatyi-Kogan & Lovelace (1997). In this ease, the 
temperature difference is expected to remain rather small 
because the heating is effective for effective radiator. In 
any case, the resistive ADAF model in its present version 
assume a common temperature to both components, for 
simplicity. The examination of its two-temperature ver- 
sion may belong to a future work. 



3. CALCULATION OF SPECTRUM 

As mentioned above, the radiation spectrum from a re- 
sistive ADAF is calculated based on the analytic solution 
introduced in the previous section. Back reactions of the 
radiation cooling to this fully advective solution are negli- 
gible as far as its fraction in the total cooling rate is small. 



This has been roughly checked in a previous paper (KOO) . 
The method of calculating radiation fluxes described in 
this section will be applied to Sgr A* in the next section. 
The observed spectrum of Sgr A* in the frequency range 
from radio up to X-ray range is successfully explained in 
the viscous ADAF models by the three processes, i.e., syn- 
chrotron radiation, bremsstrahlung and inverse Compton 
scattering (Narayan, Yi & Mahadevan 1995; Manmoto et 
al. 1997; Narayan et al. 1998). Although there may be 
some other components such as the radio-frequency excess 
over the Rayleigh- Jeans spectrum and 7-ray peak both of 
which need separate explanations (see, e.g., Mahadevan 
1999 for the former, and Mahadevan, Narayan & Krolik 
1997 for the latter), we ignore these components here for 
simplicity. 

Among the above three processes, the Compton scat- 
tering is treated separately from the other processes of 
emission and absorption. Therefore, we divide the total 
flux into two parts: the flux due to bremsstrahlung and 
synchrotron process and that due to the inverse Comp- 
ton process F^^. The obtained fluxes are both integrated 
over the entire surfaces (upper and lower ones) of a disk, 
and added up to obtain the luminosity per unit frequency 

Temperature in the flow is vertically isothermal in the 
present model. In calculating the emission and absorption 
processes, the flow is further assumed to be locally plane 
parallel. Solving the radiative diffusion equation at a given 
radius i?, we obtain the flux of the unscattered photons Fi, 
emanating from one side of the disk (Rybicki & Lightman 
1979) as 



2tt 



1 - exp(-2\/3T;) 



(10) 



where B^, is the Planck intensity and r* is the vertical 
optical depth for absorption. 



-Ki,RA. 



(11) 



Assuming the local thermodynamic equilibrium (LTE), 
we can express the absorption coefhcient k^, at the equa- 
torial plane in terms of the volume emissivities Xi^'s for 
bremsstrahlung and synchrotron processes: 



(12) 



3.1. Bremsstrahlung 

At relativistic temperatures, we must take into account 
not only electron-proton but also electron-electron encoun- 
ters. Therefore, the total bremsstrahlung cooling rate per 
unit volume is written as 



9br = 1e 



(14) 



where the subscripts ei and ee denote the electron-ion 
and electron-electron processes, respectively. The explicit 
expressions of the cooling rates are as follows. For the 
electron-proton process, 

= 1.25 nlaTcafrrieC^ Fei{9e) 

= 1.48 X 10-^^nlFe,{ee) ergs cm'^ s'^ (15) 

where n,, is the electron number density, a/ is the fine- 
structure constant and ctt is the Thomson cross-section, 
and further 



0.5 



(1 + 1.7816i,'-^'*) 



for 



< 1, 



9^ 

27T 



[ln(1.1236ie + 0.48)4- 1.5] 



for de > 1. 
For the electron-electron process, 
20 



(16) 



Qee = nlcrlm^c^af 



(44 - 37r2)6l^/2 



x(l + 1.16le + 
= 2.56 X IQ-'^'^n'^ 



97rO-5 

el - 1.2561^/2) 



6*3/2(1 + I.l6le + 6|2- 1.256*^/2) 



ergs cm ^ s ^ 



(17) 



when 6'p < 1, and 



= nlcrlmeC^af2We{\n2'qee + 1.28) 
= 3.40 X 10-22n26'e(lnl.l236le + 1.28) 

ergs cm^"^ s^^ 



(18) 



when 6*2 > 1. Here, — /rUeC^ is the classical electron 
radius and rj — exp(— 7^) = 0.5616. 

The emissivity per frequency is given by 



Thus, equation (M) includes not only the effect of free-free 
absorption but also of synchrotron self-absorption at low 
frequencies. 

As the distribution function for thermal electrons, we 
assume that of the relativistic Maxwellian (in its normal- 
ized form). 



7^/3exp(-^) 

()eK2{i:) 



d7. 



(13) 



because ADAFs tend to have so high temperatures that 
electron thermal energy can exceed its rest mass energy. 
Here, 7 is the Lorentz factor, fee is the Boltzmann constant 
and K2 is the 2nd modified Bessel function. Actually, we 
use this formula only in the calculation of Comptonized 
photon flux below, while in those of bremsstrahlung and 
synchrotron processes we follow the works of Narayan & Yi 
(1995b), and Manmoto et al. (1997) where it is replaced 
by a numerical fitting function. 



XuM = gbr^ exp 



hv 



ergs cm 



-iHz- 



(19) 



where h is the Planck constant and G is the Gaunt factor 
which is written (Rybicki & Lightman 1979) as 



G 



k^Te hv 
k^Te TT \C, hv 



for 
for 



hv 

k^Te 

hv 
kuTe 



> 1, 
< 1. 



(20) 



The above cited formule contain a few minor de- 
fects. For example, the non-relativistic limit calculated 
for electron-ion process from equations (|l^) and ( |20| ) dif- 
fers by about 35% from the standard formula (RyEicki & 
Lightman 1979). Equation (EOh assumes the same values 
of the Gaunt factor for both electron-electron and electron- 
ion processes. In spite of these defects, we adopt the above 
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formule according to Narayan & Yi (1995b) and Manmoto 
et al. (1997), considering that these are the best ones 
we can employ at present throughout the energy rage of 
our interest. The adoption of the same formule as in the 
previous calculations is also suitable for the purpose of 
comparison of the predictions of different models, such as 
the viscous and resistive ones. 

3.2. Synchrotron Emission 

Synchrotron emission is an essential process to produce 
the radio wave-length part of the spectra from optically 
thin ADAFs in AGNs. Especially in the resistive ADAF 
model, some information about the strength of the ambi- 
ent magnetic field may be obtained from the process of 
spectral fitting. 

The optically-thin synchrotron emissivity by relativis- 
tic Maxwellian electrons is calculated from the formula 
(Narayan & Yi 1995b; Mahadevan, Narayan & Yi 1996), 



X.,sy = 4.43 X 10 



_3Q A-KrieV I ( AnmeCi' 

K2{i/e,) [ ieBei 

-3 „-l 



ergs cm " s " Hz ^ 
where e is the elementary charge and 



(21) 



, 4.0505 / 0.4 0.5316\ , „„„„ w.^^, 
^ (^) = V + exp(-1.8899xi/?^2) 



In equation (|2l|), the argument of /' is specified as 



2iy 
3^^' 



I^Q 



e\B\ ^ 

2'KmeC 



(23) 



where B is the local value of magnetic field for which we 
substitute b^. 

3.3. Inverse Compton Scattering 

The soft photons whose flux is given by equation ( p^ are 
Compton scattered by the relativistic electrons in the flow. 
We adopt the rate equation of Coppi & Blandford (1990) 
as the basis of our considerations. This equation applies 
to homogeneous, isotropic distributions. The first term on 
the right-hand side of their equation describes the rate of 
decrease in the photon's number density with a given en- 
ergy owing to the scattering into other energies, while the 
second term does the increase owing to the scattering into 
this energy from other energies. 

In the situations of our interest, we can neglect the first 
term because the number density of Comptonized photons 
are small compared with that of the seed photons. Instead, 
we use the second term iteratively to calculate the effects 
of multiple scattering. The scattering occurs on the av- 
erage when the condition caTUedt = 1 is satisfied, where 
t is time and is the number density of electrons. The 
probability that such a condition is satisfied j-times before 
the photons come out of the surface may be given by the 
Poisson formula. 



e 



Pj = 



(24) 



Then, the production rate for the photons with a normal- 
ized energy e = hv/mf.c^ is given by 



dn{t) 

CaTUedt 



oo ^ 

Y,Pj / djj ...d-fi iV,(7,) . . . iVe(7i) 
i=i 



dcj . . . dei 



P(e; ej,jj) . . . P{e2] ei,7i) i?(ej,7j) . . 



i?(ei,7i) nin(ei) 



(25) 

is the number den- 



where me is the electron mass and n-^i 
sity of seed photons. 

The non-dimensional scattering rate R{e,"f) includ- 
ing Klein-Nishina cross section ctkn is written explicitly 
(Coppi & Blandford 1990) as 



(TKNiP, e, m) 



27(l+/3)£ 



3272/3e2 72^(,_^)^ 

1 8 1 
^2^x^ 2(H-x)2 



dx 



^]Hi + x) 

X X 



s — 1 

cm s . 



(26) 



Scattered-photon distribution is denoted by P(e; e', 7) and, 
in the present calculation, approximated by a 5-function 
(Lightman & Zdziarski 1987, Fabian et al. 1986): 



(27) 



This is merely for simplicity and a more exact expression 
has been derived by Jones (1968) and corrected afterwards 
by Coppi & Blandford (1990). 

Although the sum in equation (|25| ) runs to infinity, it 
seems appropriate to assume that photons which are scat- 
tered more than certain times become saturated and obey 
the Wien distribution cx exp{— hiy / kT^) (e.g., Manmoto 
et al. 1997). In view of the smallness of the optical depth 
in most sub-Eddington ADAFs (res < 10~^, in the case of 
Sgr A*), however, we truncate the power series in Tes at 
j — 2 and ignore the saturation effect. After performing 
the integrations containing (5-functions and transforming 
the photon number densities into fluxes by multiplying 
che/2 on both sides of equation (Ea), we obtain 



-F. 



(28) 



where once- and twice-scattered fluxes are given, respec- 
tively, by 



d7i 7Ve(7i) R 



3e 



71 



3e 



^72 



1 

X R 

X F,- 



47i / V47i 
d-ii Nei-/2) iVe(7i) 



(29) 



3e 



R 



3e 



472^ 47i^ 



71 



3e 



472' 47/ 



(30) 



with the definition f^jn = (c/ie/2)nin(e) = Jin(e)- As the 
incident flux F^^in in the above expressions, the result from 
equation (nW should be used. The effects of low energy 
tail (/3 < 1) in the electron distribution is neglected in 
performing the 7-integrals. 
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4. APPLICATION TO SAGITTARIUS A* 

In order to compare the resistive ADAF model with the 
current models of viscous ADAFs in their predictions of 
spectra from accretion flows, we apply the former model 
to Sgr A* . The observed spectral data available so far have 
been compiled by Narayan et al. (1998). They assume that 
the interstellar column density is A^h = 6 x 10^^ cm~^ 
and the distance to the Galactic center is d = 8.5 kpc. 
In judging the accuracy of fittings between the calculated 
and observed spectra, a considerable weight has been put 
on the high resolution data points such as the VLBI radio 
(86 GHz) data and the ROSAT X-ray data. It should be 
kept in mind, however, that the ROSAT data may be in- 
terpreted as an upper limit because its resolution (PSPC) 
~ 20" is not considered as satisfactory and that other is- 
sues like the value of A^h are still under discussion. 

We discuss the two cases in the resistive model, which 
will be called the compact-disk and the extended-disk 
models, respectively. These names come from the differ- 
ence in extension of the disk which is represented by the 
radius ratio of the inner to the outer edges, -Rin/^?out = Hn- 
As confirmed below, this value is largely affected by the 
choice of position of the inner edge, Rin- 

4.1. Com,pa,ct Disk Model 

According to the spirit of original resistive ADAF model, 
the inner edge in this case is determined by the magnetic 
flux conservation (KOO). This gives an expression 



Ri, 



(l + A-i)-2i?,,t~ A^i?, 



(31) 



where the last expression is valid only for thin disks 
(A <C 1). Note that this procedure is independent of the 

notion of the marginally stable orbit around black holes. 
The outer edge has been fixed, on the other hand, from 
the mass conservation as 



Ro 




(32) 



Fig. 2 shows the best fit spectrum in this model and the 

set of best fit parameters is 



M = 3.9 X 10'^ Mq, 

M = 1.2 X IO^^'Me = 1.0 X IQ-^Mq yr" 

|Bo| = 0.7 G, A = 0.14 rad. 



(33) 



From these values, other quantities of our interest are fixed 

as follows: 



Rin = 6.1 X 10^^ cm, 
T = 3.4 X 10^ r"^ K, 



^ 3.1 X 10'* cm, 
\b^\ = 5.0 r-^ G, 
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p= 1.8 X 10 

Tes = 3.1 X 10"^. 



g cm 



(34) 

Thus, it turns out that the inner edge of the present model 
is fairly large compared with the marginally stable orbit, 

i?nis = 3.5 X 10^^ cm, for a Schwarzschild hole of the above 

mass. 

The changes in the spectrum caused by varying cen- 
tral mass M, accretion rate M, external magnetic field 
\Bq\ and disk's half-opening angle A are demonstrated in 
Figs. 3, 4, 5 and 6, respectively. The spectral features 
are anyway quite analogous to those predicted by the vis- 
cous ADAF models. The results of a detailed comparison 
between the resistive and viscous ADAF models will be 
discussed in the final subsection, based on the predicted 
spectral features. 



4.2. Extended Disk Model 

In this model, the inner edge of the accretion disk is set 
at the radius of the marginally stable circular orbit around 
a Schwarzschild black hole, 



Rh 



R„ 



3i?G 



6GM 



(35) 



where Rq is the gravitational radius of the hole. This 
choice is motivated by the expectation that at around this 
radius the infall velocity inevitably becomes of the order 
of the rotational velocity, (i.e., 3? ^ 1 where 3? is the mag- 
netic Reynolds number, see K99, KOO). The definition of 
the outer edge is the same as in the compact disk model. 
Note that the above definition of the inner edge is adopted 
also in the viscous ADAF models. 

The best fit parameters in this model are 

M = 1.0 X 10^ Mq, 

M = 1.3 X IQ-^Me = 2.9 X 10"^ Mq yr'^ 

\Bo\ = 1.0 X 10"*^ G, A = 0.20 rad. (36) 

These are used to fix the values of various scaled quanti- 
ties: 



Rin = 8.9 X lO" cm, 
r=1.0xl0*r"^ K, 
/> = 5.9 X 10"^^ r"^ g cm" 

Tes = 1.3 X 10"^ 



Rout = 2.7 X 10 
16, 
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5.0 X 10"*^ r"^ 



G, 



(37) 



The best fit curve is shown in Fig. 7. The changes in 
the spectrum caused by varying central mass M, accre- 
tion rate M, external magnetic field |i?o| and disk's half- 
opening angle A are demonstrated in Figs. 8, 9, 10 and 11, 
respectively. The spectral shapes are very different from 
those of the compact-disk case and of the viscous ADAF 
models. Synchrotron emission has a very wide peak and 
bremsstrahlung is negligibly small. The former fact is due 
to a high temperature at the inner edge (see sub-subsection 
4.3.2) and the latter, to lower densities in the disk. The 
emission in the X-ray band is supported by the inverse 
Compton scattering from the radio band. The tempera- 
ture near the outer edge falls even to such a small value 
that the assumption of complete ionization becomes in- 
valid. Although the position of outer edge may seem to be 
irrelevant from a viewpoint of spectrum, it is nevertheless 
important also in this case as a fitting boundary of the in- 
ner magnetic field to the external one. The fitting predicts 
that the boundary value is comparable to the interstellar 
field (a few /iG). 

The fitting both to 86 GHz and ROSAT X-ray data 
points is possible also in this model. However, it is clear 
that the fitting curve runs above the observed upper limits 
in the IR band. The fitting in the frequency range from 
100 to 1000 GHz also becomes considerably poor com- 
pared with the case of compact disk. For these reasons, 
we judge that this model cannot reproduce the observed 
broadband spectrum of Sgr A*. This fact suggests again 
that the inner edge of the accretion disk does not coincide 
with the marginally stable orbit. The wide range of the 
disk's radii which is obtained from this fitting implies that 
^{Rout) ~ 6 X 10^. Since ^{R) represents the ratio of 
toroidal to poloidal magnetic fields, most parts of the disk 
are very likely to be unstable to global MHD instabilities 
of helical type. For this reason too, we consider that the 
present case (i.e., Rin = Rms) is quite unrealistic, at least, 
for Sgr A*. 
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4.3. Viscous v.s. Resistive ADAFs 
4.3.1. Dependence on Black Hole Mass 

The spectra calculated from ADAF models of both vis- 
cous and resistive types commonly have the saturated part 
at the lower ends of the spectra due to the synchrotron self- 
absorption. It is of great interest to see that the luminosity 
vLi, of this part is essential to determine the mass of the 
central black hole, in both types of models. Especially, in 
the viscous model, the luminosity of this frequency part 
is determined almost only by the black hole mass. The 
reason is as follows. 

The temperature in ADAFs may be considered essen- 
tially as the ion virial temperature and hence decreases 
as ~ Apart from a numerical factor due to a re- 

duced Keplerian rotation, this is exactly true in the re- 
sistive model. This is also true in the viscous models for 
the main part of an accretion flow except in the inner re- 
gion where the electron temperature deviates from the ion 
temperature and remains almost constant (e.g., Narayan 
& Yi 1995b). Therefore, the contribution to the spectrum 
from each annulus of radius R and width dR is equal. 
Integrating these contributions up to the outer edge, we 
obtain Lf^ cx Te{Rin)RinRout = Te(i?out)^out> where i?in 
is the radius of the disk's inner edge in the resistive model 
and of the outer edge of the two-temperature region in the 
viscous models. 

We have T^R (x m commonly to both types of ADAF 
models. Further, since radius scales as the gravitational 
radius in the case of viscous ADAFs, we obtain the mass 
dependence 



(viscous ADAF), 



(38) 



confirming the above statement. On the other hand, in 
the case of resistive ADAFs, we have 



(resistive ADAF), (39) 



where the dependences on the parameters other than m 
have come from the expression of i?out- In spite of these 
dependences, the mass dependence is essential also in this 
case. This is because the dependence on m is rather weak 
and the value of bo is strongly restricted from the position 
of the synchrotron peak (see the discussion below). 

4.3.2. Synchrotron Peak 

We estimate the synchrotron peak frequency following 
Mahadevan (1997), and examine its behavior in both vis- 
cous and resistive models. For each annulus of radius R 
and width dR, the synchrotron photons in the radio range 

up to a critical frequency are strongly self-absorbed and 
result in the Rayleigh- Jeans spectrum. Therefore, the crit- 
ical frequency of the spectrum is determined by equating 
the contributions to from optically thick and thin sides 
of the frequency: 



2Tr-^kBTe(R)2TTR dR 

= 4.43 X 10-30 -^^^/'(a;e) ^nAR^ dR, (40) 

A2(l/c'e) 

where Xc is defined as Xc = {'i^oOe)- Solving this equa- 
tion, we can determine the value of Xc numerically (Ap- 
pendix B of Mahadevan 1997). Provided that this value 
does not depend strongly on R, A and other parameters, 
we obtain 



If the disk has uniform temperature and magnetic field, 
then the synchrotron peak is rather sharp and has a well- 
defined peak frequency at Vc- 

When they vary with the radius R, however, substitu- 
tion of the r-dependences of To and B in both viscous and 
resistive ADAF models yield 



V, cx _ /3)l/2^1/2^-l/2^-13/4 

(viscous ADAF), 



(resistive ADAF) . 



(42) 



This means that is larger for smaller radii and the 
higher most cutoff is due to the inner edge. The posi- 
tion of peak of the superposed emission is then given as 
t^p = i^c{f-p), where rp is the radius whose contribution 
to the synchrotron emission is most dominant. The fairly 
narrow peak obtained in the compact-disk case indicates 
that Tp is located near the inner edge and the global peak 
shape is determined mainly by the inner most region of 
the disk. 

On the contrary, the synchrotron peak becomes very 
broad and dull in the extended-disk case. We have con- 
firmed that the low-frequency side of the broad peak is 
due to a superposition of the contributions from annuli 
of i?ms ^ 10i?ms- However, the dull shape on the high- 
frequency side of the peak may be mainly due to a re- 
sulting high temperature (T ~ 3 x 10^" K) at the smaller 
inner edge. Actually, owing to this high temperature and 
low densities near the inner edge, the synchrotron self- 
absorption becomes less important in the high-frequency 
radio band and the intrinsic shape of the synchrotron emis- 
sion at the mildly relativistic temperature (Mahadevan et 
al. 1996) can appear on the high-frequency side. 

In any case, since rp is a numerical factor, we can speak 
of the parameter dependences of the peak frequency 



based on equation (42). Note that the dependences on 

m and m have different senses in the different ADAFs. 
The most important difference between the two models is 
that the dependence of Vp on the magnetic field is much 
more sensitive in the resistive model. Therefore, the field 
strength is determined more accurately there. All the pre- 
dicted dependences on to, m, 5o and 5 are qualitatively 
confirmed in Figs. 2 through 5. From the above consid- 
erations on the synchrotron peak, we think that the im- 
provement of observational quality in submillimeter range 
is most important for obtaining more exact values of the 
disk parameters. 



4.3.3. Bremsstrahlung 

We shall try here to grasp the qualitative behav- 
ior of the contribution from bremsstrahlung according 
to the usual non-relativistic scheme.^ The , contribu- 
tion to a given frequency v from optically thin plasma 
in an annular volume of width dR is proportional to 
p^T~^/'^ eyip[—hv /kBT] R^dR. Apart from the exponen- 
tial factor, we have 



o^T ^^^R^dR cx a m mr 



= flv^x, (X T^\r)B{r). 



'^rri'mr ^^^dr (viscous ADAF), 

jc-4i,-2/5 . 11/5 4/5 1/2 1 

CX Oq To ' To ' r ' dr 



(41) 



(resistive ADAF). 



(43) 
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Therefore, the relative importance of the inner and outer 
parts of a disk can be seen from the ratio, 



(44) 

whore C = -Rout / -Rin = J'in^ , and the upper and lower signs 
in its exponent are for the viscous and resistive ADAFs, 
respectively. 

It is evident from the above ratio that, in the viscous 
ADAF, the contribution from the inner disk is always dom- 
inant (i.e., / > 1) irrespective of the frequency v. In the 
resistive ADAF, however, it depends on the frequency, so 
that we introduce the critical frequency v'^ by the relation 
/ = 1. This yields 



2 h 



(45) 



Then, the inner part contributes to the frequency range 
V > v'c and the outer part, to v < v'^. In fact, the critical 
frequency roughly coincides with the peak frequency of the 
bremsstrahlung. The luminosity above can be roughly 
estimated, by putting dr r 

7*in, clS 

L^'' oc a~^m^m (viscous ADAF), 

oc (resistive ADAF), (46) 

because rin is a numerical constant and, in particular, 
equal to S'^ in the resistive ADAFs. 

In the viscous ADAF models and in the compact- 
disk case of the resistive ADAF, the contributions from 
bremsstrahlung cause an X-ray bump in each predicted 
spectrum. The dependences of this peak on the parame- 
ters m and rh are qualitatively confirmed in Figs. 2 and 
3, but those on bo and 6 are somewhat different from the 
above prediction, indicating a limitation of such a crude 
estimate as the above. The critical frequencies calculated 
from the best fit values for the compact and extended disks 
arc 1.4 X lO^*^ Hz and 1.8 x 10^'"^ Hz, respectively. The 
former value is in good agreement with the peak of the 
reproduced spectrum. In the extended-disk case, the con- 
tribution from the bremsstrahlung is negligibly small be- 
cause the density throughout the disk becomes too small, 
and the X-ray range of the spectrum is explained by the 
once-scattered Compton photons. 

5. SUMMARY & DISCUSSION 

To summarize the examinations in the previous section, 
both viscous and resistive ADAF models can explain the 
observed spectrum of Sgr A* equally well. In spite of large 
differences in the basic mechanisms working in both mod- 
els, the calculated spectra are quite similar, except for the 
extended-disk case in the resistive model. This fact sug- 
gests that also the resistive ADAF model is quite powerful 
in explaining the behavior of other low luminosity AGNs 
(Narayan, Mahadevan & Quataert 1998). In addition to 
these analogous aspects, the resistive model seems to have 
a possibility to explain such an esscnitially different situa- 
tion as appeared in the extended disk case. In any case, 
when the presence of an ordered magnetic field should be 
taken seriously in some AGNs or in some stellar-size black 
holes then the resistive ADAF model, whose predictions 
on the radiation spectra are examined in this paper, will 
serve the purpose. 

One of the most remarkable features of the ADAF mod- 
els is that the mass of the central black hole seems to be 



determined only from the fitting to the self-absorbed part 
of the observed spectrum. In the case of Sgr A*, the re- 
sistive ADAF model (hereafter restricting to the case of 
compact disk) predicts the central mass of 3.9 x 10^ Mq 
while the viscous models predict 1.0 x 10^ Mq (Manmoto 
et.al 1997) and 2.5 x 10^ Mq (Narayan et al. 1998). The 
accuracy of the fittings for other disk parameters than the 
black hole mass will be greatly improved by the precise de- 
termination of the position and height of the synchrotron 
peak from observations. 

The black hole mass predicted by the resistive ADAF 
model is evidently smaller compared with the predictions 
of the viscous ADAF models. The latter values are consis- 
tent with the dynamically reduced value of 2.5 x 10^ Mq 
(Haller et al. 1996; Eckart & Genzel 1997), which may be 
considered as an upper limit for the black hole mass. In 
the history of viscous ADAF models, the predicted black- 
hole mass was as small as 7 x 10^ Mq (Narayan et al. 
1995). Afterwards, by the inclusion of compressive heat- 
ing, it becomes consistent with the dynamical mass. Since 
this change is mainly due to the decrease in electron tem- 
perature (Narayan et al. 1998), the prediction of the resis- 
tive ADAF model may also bo increased if the development 
of its two-temperature versions results in a lower electron 
temperature. As for the compressive heating, it is already 
included in the resistive model. 

In spite of the resemblance in the predicted spectral 
shape, there are of course many differences in the pre- 
dictions of the viscous and resistive models. The precise 
dependences on the relevant quantities of the luminosities 
of the self-absorbed part, the synchrotron peak and the 
X-ray bump arc different. Especially, the dependence of 
the synchrotron peak-frequency on the strength of mag- 
netic field is much stronger for the resistive model. The 
essential difference in the geometry of an accretion flow 
may be in the radius of the inner edge rather than in its 
vertical thickness. The prediction of the resistive ADAF 
model for the inner-edge radius of the disk around Sgr 
A* is ^ 20i?ms5 instead of the radius of marginally stable 
circular orbit i?ms- Although this result justifies the ne- 
glection of the general relativistic effects in our treatment, 
various questions may be raised about the behavior of in- 
falling plasmas. As for this point we only present an idea 
below from a viewpoint of global consistency, because its 
detailed analyses arc beyond the scope of this paper. 

Fig. 1 shows an overview of the flow and magnetic field 
configurations (see KOO, for more details). The accretion 
flow would be decelerated near the inner edge by the pres- 
ence of a strong poloidal magnetic field which is main- 
tained by the sweeping effect of the flow. As a result, a 
certain fraction of the accreting plasma will be turned its 
direction to go along the poloidal field lines, although the 
remaining fraction may fall into the central black hole. If 
the poloidal current driven in the accretion disk can close 
its circuit successfully around distant regions and along 
the polar axis, a set of bipolar jets will be formed (Kabu- 
raki & Itoh 1987). Even if the mechanism for formation of 
jet does not work well, the plasma within the inner edge 
is likely to extend to the polar regions. 

The presence of the plasma within the inner edge of an 
accretion disk and near the polar axis can be a possible 
source of the excess above the self-absorbed slope in ra- 
dio band of the observed spectrum. Very recent VLBI 
observations of Sgr A* (Krichbaum et al. 1998; Lo et al. 
1998) report that its intrinsic sizes in the east-west di- 
rection at 215 GHz and 68GHz are about 20 i?G (with 
M = 2.5 X 10'^ Mq). a haff of this size (i.e., its radius) 
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is just comparable to the size of the inner edge ~ 60i?G 
of our model fitting with M = 4 x IQ'^Mq. However, the 
value of the black hole mass estimated from the spectral 
fitting may be increased if there is a possibility for Sgr 
A* to have a wind-type mass loss from the surfaces of its 
disk (such possibilities have been noted for various types 
of objects by, e.g., Blandford & Begelman 1999; Di Matteo 
et al. 1999; Quataert & Narayan 1999). In such a case, 
the VLBI component becomes smaller than the size of the 
inner edge. 

From the standpoint of the resistive ADAF, therefore, 
the above observations should be interpreted as suggesting 
the presence of a compact structure which is comparable 
to or smaller than the inner-edge radius of the accretion 



disk. In view of the vertical elongation of this component 
reported by Lo et al. (1998), this structure is very likely 
to be the root of a jet as suggested by them. This pic- 
ture is very consistent with the view described above in 
relation to Fig. 1. In this case, however, the location of 
the self-absorbed slope in the v-vL^ diagram should be 
slightly shifted towards the higher-frequency side so that 
the VLBI data points can be regarded as an excess from 
the disk's contribution. 

One of the authors (M. K.) would like to thank Tadahiro 
Manmoto for many valuable comments on the viscous 
ADAF models. He is also grateful to Umin Lee for his 
suggestions on some numerical technics. 
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return current 



Fig. 1. — Schematic drawing of the global geometries of magnetic field and plasma flow. A poloidally circulating current system (jp) 
driven by the rotational motion of accreting plasma generates a toroidal magnetic field b,p in addition to a nearly uniform external field. The 
presence of this toroidal field outside the disk guarantees the magnetic extraction of angular momentum from the disk. This field also acts to 
confine the accreting flow toward the equatorial plane and has a tendency to coUimate and accelerate the plasma in the polar regions. If the 
condition is favorable, the plasma in the polar regions may form a set of bipolaj: jets. In this paper, however, we focus our attention on the 
radiation spectrum from the accretion disk only. 
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Fig. 2. — The best fit spectrum of Sgr A* in the compact disk model. The resulting physical quantities arc M = 3.9 X 10^ A^q, 
M = 1.2 X 10^* A^Ei l^^ol = 0.7 G and A = 0.14 rad. The four peaks indicated by S, CI, C2 and B denote those due to synchrotron emission, 
once- and twice-scattered (although it is almost buried) Compton photons and bremsstrahlung, respectively. The data points are the same 
as those compiled by Narayan et al. (1998). 
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Fig. 3. — With increasing central mass, the luminosity vLi, increases globally, so that the self-absorbed part shifts upward and the 
synchrotron pealj-frequency sifts to higher frequencies. The direction of the shift of the synchrotron peak is different from the results of 
viscous ADAFs (Manmoto et al. 1997, Narayan et al. 1998). 
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Fig. 4. — M-depcndcncc of the spectrum. With increasing mass accretion rate, tiie synchrotron peak shifts to lower frequencies while the 
self-absorbed luminosity is hardly affected. The former behavior is opposite to that of viscous ADAFs. The X-ray bump increases with M. 
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Fig. 5. — |Bo|-dependence of the spectrum. With increasing magnetic field strength, the synchrotron peak shifts to higher frequencies while 
the self-absorbed luminosity is hardly affected. The X-ray bump increases with increasing |Bo|; suggesting a limitation of the crude estimate 
of iJ^ given in the text. 
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Fig. 6. — A-dependence of the spectrum. With increasing half-opening angle, the synchrotron peak shifts to lower frequencies while the 
self-absorbed luminosity is hardly affected. The dependence of the X-ray bump on A again suggests a limitation of the crude estimate given 
in the text. 
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Fig. 7. — The best fit spectrum of Sgr A* in the extended disk model. The resulting physical quantities axe M = 1.0 X 10^ Mq, 
M = 1.3 X 10"'* Me, \Bo\ = 1.0 x 10~® G, A = 0.20 rad. Among the four peaks which are seen in the case of compact disk, that of 
bremsstrahlung has disappeared owing to the low densities in the extended disk, so that the X-ray bump is explained by once-scattered 
Compton photons. The fitting both to 86 GHz and ROSAT X-ray data points is possible also in this model. However, the fitting in the 
frequency range from 100 to 1000 GHz becomes considerably poor compared with the case of compact disk, and it cannot be reconciled with 
the observed upper limits in the IR band. For these reasons, we judge that this model cannot reproduce the observed broadband spectrum 
of Sgr A*. 
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Fig. 8. — M-dependence of the spectrum. The tendency of the changes caused by varying M is almost the same as in the case of compact 
disk. 
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Fig. 9. — M-dcpcndcncc of the spectrum. With increasing mass accretion rate, the synchrotron peak shifts to higher frequencies while the 
self-absorbed luminosity is hardly affected. The former tendency is different from the case of compact disk. 
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Fig. 10. — |Bo|-dependence of the spectrum. The tendency is almost the same as in the case of compact disk. 
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Fig. 11. — A-dependence of the spectrum. The tendency is almost the same as in the case of compact disk. 



